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consistent mass approach over estimates frequencies while lumped 
mass approach under estimates. Difference could also be due to the 
fact that bearing properties used in the present work may be 
different from those of C17D (not given). 
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properties as shown in Table 13. 
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From the unbalance response of shafts, the following 
observations are made: 

(a) With no damping the system becomes unstable at critical 
speed . 

(b) External damping when introduced, either from medium or from 
bearing, the system stability is increased and whirl amplitudes 


decreases . 
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ABSTRACT 

Flexural vibrations of damped shafts including the effect of 
shear deformation, rotary inertia and gyroscopic action are best 
represented by four second order differential equations. Transfer 
finite element method (TFEM) has been developed in the present 
work for determining its critical speeds, complex natural 
frequencies and unbalance responses. By using four second order 
differential equations, all types of geometrical and natural 
boundary conditions, including bearing cross coupled properties 
can easily be incorporated exactly. This method blends the 
advantages of standard finite element method and transfer matrix 
method, i.e* solution over an element is obtained in a straight 
forward manner and memory requirements are considerably reduced. 
A variety of situations are considered and results are found to be 


satisfactory. 



CHAPTER I 


INTRODUCTION 


1.1 THE PROBLEM IN BRIEF 

Dynamics of rotating shafts, rotor dynamics in short is a 
very important area in the field of engineering as rotors are an 
integral part of many machines like steam and gas turbines, 
compressors, generators, motors etc. 

It was recognised early in eighteenth century that a rotating 
shaft initiated vibrations of large amplitudes at certain speedy 
and often caused severe damage to the machinery system. The 
operating speedSat which this phenomenon occured came to be known 
as critical speeds. 

A well balanced shaft also has a certain amount of unbalance 
which causes it to deflect from its designed position during 
rotation. When the shaft rotates in the deflected position, the 
disturbing inertia forces rotate with the shaft and shaft 
precesses about the bearing axis. This precesion of the shaft can 
be in the direction of rotation called forward precession or in 
the opposite direction called backward precession. 
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The inertia force due to unbalance has the frequency of spin 
speed of the shaft- Hence when the speed of precesion of shaft 
becomes equal to the spin speed , resonant state is reached, i.e. 
amplitude grows with time. When this state is reached, the 
corresponding spin speed of the shaft is called critical speed - 
Depending on the direction of precession if is called forward or 
backward critical speed. 

1.2 LI TERATURE RE VI EW 

The earliest work published relating to this problem was in 
1869 by Rankine Ell. The operating speeds of rotating machines 
then were much below the first critical speed. Hence most of the 
work then were limited to the evaluation of first critical speed. 
But with advent of new technology operating speeds have became 
higher than the first critical speed. Consequently more recent 
studies treat a greater range of problems with study of effect of 
various system parameters on critical speeds, determination of 
stable range of operation etc. In 1967, Loewy and Piarulli L21 
tried to illustrate the basics of rotor dynamics, in the monograph 
Dynamics of rotating shafts. This work treated shaft as mass less 
i.e. used lumped mass approach- It provided a general but 
qualitative information on the effects of various parameters on 
critical speeds and stability- This book also includes an 
exhaustive bibl iography, 

Tondl E33 in his book gives a good insight into several 
phenomenon of rotor dynamics. He emphasises on effect of internal 
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damping, unequal shaft stiffness and oil film in journal bearings. 
It deals comprehensively with problems concerned with the 
stability of rotor motion, with the creation of self excited 
vibrations and with nonlinear resonances . But almost the whole 
book is based on the assumption of mass less shaft and rigid disk. 
All the same the book gives an excellent review of the problems, 
supplemented with results of experimental investigations. 

Dimentberg £43 deals with flexural transverse vibrations in 
his book. Major part of the book uses continuous approach. 
Effects of external and internal damping, dissimi larity of 
stiffness, multi bearing systems, uniform stiffness and 
non-uniform stiffness shafts with and without shear deformation 
effects are studied in this book. Effect of shear deformation is 
studied quite in detail in this text book. 

Two other recent texts in this field are of Rao £53 and 
Goodwin £63. These books deal with various aspects of rotor 
dynamics using lumped mass approach. However, the effect of 
rotary inertia and gyroscopic action of the disks (modelled rigid) 
have been considered but neglected for shafts. Another recent 
text is Lalanne and Ferrar is £73. This book uses continuous 
approach and Rayleigh Ritz finite element method. In these texts 
£5,6,73 the effect of shear deformation is not discussed. 

An excellent collection of papers on several aspects of rotor 
dynamics, presented at an international symposium held at Denmark 


in 1974 can be found in the book edited by Niordson £83. 


Another 
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such collection of papers are found in "Vibration and wear in high 
speed machinery", edited by Silvas £91 . 

The application of finite element method (FEM ) to the 
problems relating to rotor dynamics was done by Ruhl and Booker 
£101, in 1972. They used Euler beam theory incorporating only 
translational kinetic energy and elastic bending theory to analyse 
the stability and unbalance response of a turbo rotor system. 
Nelson £113 improved the model by incorporating effects of rotary 
inertia and gyroscopic action. The effects of bearing stiffness 
and damping are also incorporated . Zorzi and Nelson £ 123 , 
incorporated the effect of internal damping in the finite element 
model of £111. The mechanism of internal damping is better 
represented by considering the hysteretic form rather than viscous 
form £41. Reference £121 studied a model combining both the forms 
of damping. They also studied effect of bearing stiffnesses and 
damping charac ter ist ics on the instability threshold. Later in 
1980, Nelson £13 1 used FEM, incorporating effect of shear 
deformation in addition to rotary inertia and gyroscopic action. 

Mahendrakar £141 studied the critical speeds, natural 
frequencies and stability of rotating shafts by Galerkin FEM, 
including the effects of rotary inertia, shear deformation and 
gyroscopic action. He used two simultaneous second order 
differential equations in comp 1 e ;< variables for undamped systems, 
and two simultaneous fourth order differential equations in real 
variables for damped systems, but with shear deformation neglected 
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in this. He did not study effect of bearing cross-coup led 
propert ies . 

Nelson and Glassgow E151 studied natural frequencies of a 
shaft supported in elastic hydrodynamic bearings. They used 
component mode synthesis using finite element analysis, but 
neglecting effect of shear deformation. 

A good number of authors have used transfer matrix method 
with lumped mass assumption. Lund E 163 studied the effect of 
fluid film journal bearings and hysteresis internal damping on the 
critical speed and instability thresh hold. However, he treated 
the shaft as a uniform beam and neglected effect of shear 
deformation, rotary inertia and gyroscopic action. Bansal and Kirk 
C 171 improved this work incorporating effect of gyroscopic 
moments, rotary inertia and shear deformation, and bearing 
characteristics . Muller's quadratic interpolation technique was 
employed to extract complex eigen values. Ramakrishnan and Frabhu 
E18I studied the unbalance response, instability and critical 
speeds of multispan rotor systems. Doughty and Vaface E191 
studied damped torsional systems. They used Newton— Raphson method 
for determining complex eigen solutions. 

Tewari E203 in his work has developed a continuous system 
approach to analyse the sub-critical phenomenon due to gravity 
effects in rotating shafts. 

In using transfer matrix method for this kind of problems, it 
is difficult to obtain transfer matrices for rotors with rotary. 
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shear and gyroscopic effects. This difficulty is overcome easily 
in a very straight forward way in transfer finite element (TFEM) , 
Gupta, Gael, Tesar and Fillo £23,24,253. This method blends the 
advantages of finite element method and transfer matrix method. 

Goel £243 obtained the critical speeds of uniform isotropic 
undamped shaft using TFEM. The effect of shear deformation, rotary 
inertia and gyroscopic moments were included in his analysis. In 
his work two second order governing differential equations solved 
using linear polynomial led to i 1 1-condi t ioning . To avoid this 
two equations were combined into a single fourth order 
differential equation and cubic polynomials used. Although it was 
not possible to satisfy all the geometric and natural boundary 
conditions exactly in this approach, the results obtained were 
quite satisfactory. 

1.3 AIM AMD SCOPE OF PRESENT WORK 

The present work aims at obtaining the critical speeds, 
whirling complex natural frequencies, its stability and unbalance 
responses of rotating shafts by developing TFEM. The governing 
differential equations are derived* including the effects of 
gyroscopic action, rata.ry inertia and shear deformation. These 
four simultaneous second order differential equations are solved 
using Galerkin FEM. This has the advantage of satisfying all 
boundary conditions, including cross— coup led bearing properties 
exactly. The FEM equations thus developed are transformed to 
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transfer finite element matrices, and complex roots obtained by 
Newton— Raphson method . 

Studies have been done on complex natural frequencies of 
whirl of rotating shafts in rigid and elastic bearing supports. 

Hysteresis as well as viscous internal dampings have been 

% 

incorporated . 

Unbalance responses of a shaft-disk system are studied for 
different end support conditions. 



CHAPTER II 


FORMULATION OF GOVERNING EQUATIONS 

This chapter gives the governing differential equations in 
non-dimensional ised form- These governing differential equations 
have been derived by til, 14, 251, but are reproduced here for 
completeness sake . 

2. 1 SHAFT WITH NO INTERNAL DAMPING 

Governing differential equations for transverse vibrations of 
spinning shaft, incorporating the effects of shear deformation, 
gyroscopic moment and rotary inertia are derived from Hamilton's 
principle - 

Consider a differential element of length ds, as shown in 
Fig- (2.1). The kinetic energy of the shaft becomes til, 14, 253, 
as shown in Appendix. 

? i 

Te = J pA (LJ 2 + V 2 ) ds + J \ pi + S 2 ) ds 

o o — X y 

l 1 l 

+ fi-ia i ds+fioeeds 

J 2 p o ^ poxy 

o K o 
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Here pA is the mass per unit length, U and V are the translations 
along ox and oy axes, pi is the mass moment of inertia per unit 
length about the neutral axis, B^ and B^ are the bending slopes in 
ys and xs planes; I is the polar mass moment of inertia per unit 
length and is the shaft spin speed, assumed constant. Dot 
represents derivatives with respect to time. In the above 
expression <2.1.1) first term is the energy due to translation, 
second term is the energy due to rotary inertia, third term is the 
energy due to spinning of the shaft and the last term gives the 
energy due to gyroscopic moments. 

The strain energy of the shaft is 

l t i 

U = F 4- + ds + f 4 k'Gfi (ft 2 + ft 2 ) ds (2.1.2) 

e J 2 x y J 2 x y 

o J o J 

where El is the flexural rigidity; k' is the transverse shear form 

factor, a constant depending upon the shape of the cross-section , 

6 is the shear modus, A is the area of cross section, and 

are the angles of distortion due to shear deformation about y and 

x axis respecti vely . Strokes (') represent the derivatives with 

respect to s. First term represents strain energy due to bending 

and second term strain energy due to shear deformation. 

It is seen from Fig. (2.1) that the total slopes are 
IT = 6 + ft 

y y 

V' = -S + ft <2.1.3) 

X X 
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Combining eqns. (2.1.2) and (2.1.3), expression for strain energy 

U becomes, 
e ’ 

u = f % ei (e 2 + e 2 y ds 

e J 2 x y 

o 


f i k'GA [(U' - 0 > 

J 2 [ y 


2 + <v' + e^) 2 1 ds 


(2.1.4) 


The variational work done <5W due to external viscous damping 


is 


L l 

<5W = J c U <5U ds + J c V c5V ds 

o a 


(2.1.5) 


By Hamilton's principle, 

*2 

6 f' (Te - U -W) dt = 0 

t 6 

r l 

i . e . 


7 ? 

J t J {I 

t, o v 


i- I O 2 

2 p o 
| k'6A 


■ 2 ) 

4- 

1 pI 

*2 *2 
(e + e ) 
x y 



6 

e - 

1 ' 2 
4 EI <0 ~ 

' 7 

+ 0 ) ds 

“o 

X 

y 

2 X 

y 

e 

y 

) 2 

+ (V' 

* V 2 ]} 

ds dt 


'2 l 


-J J (c U 6U + c V 6V> ds dt = 0 


( 2 . 1 . 6 ) 


*1 ° 


Applying the & operator and rearranging, one gets 
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'2 l 


[- 


pAU + k'GA (U'-9 ) 

y 




6U 


pi 6 + I O 6 + (EI9')' + k'GA (U'-e ) 66 

y pox y yjy 


+ £ -pAV + k'GA (V'+e^)' - cvj <5V 


[- 


pi 6—1 Cl 6 +(EI&' ) ' — k ' GA (V'+6 > <56 Us dt 

xpoy x x | x * 


] } 


- El 6' 66 I 1 - El 6‘ 66 l L - k'GA(U'-6 >t5U| L 
x x 1 o yy*o y l o 


- k'GA(V'+6 , )<5V|^ = 0 


(2.1.7) 


The displacements <5U, <56^, <5V, 66^ are arbitrary, therefore 

pAU - k'GA (U'-S )' + ell = 0 

y 


pi 6 -1 Cl e -CEI6' )'-k'GA (U'-6 ) = 0 

y p ox y y 


( 2 . 1 . 8 ) 


pAV - k'GA (V'+S^)' + cV = 0 


pi 6 +1 Cl 6 ~(EI6' ) '+k'GA (V'+9 ) = O 

xpoy y x 


(2.1.9) 


El 6‘ 66 I 1 = O , k'GA (Li'— 6 ><5Uj l = O 
x x 'o y ‘o 


( 2 . 1 . 10 ) 


El 6' 66 I 

y y'o 


k'GA(V'+6 )<5V| = O (2.1.11) 


Equations (2.1.8) are governing differential equations in 
plane, eqns. (2.1.9) in ys plane and eqns. (2.1.10) and (2.1.1 


xs 

1) 


are corresponding boundary conditions. 
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2.2 SHAFTS WITH INTERNAL DAMPING 

References L12, 14, 163 have discussed cases of internal 
viscous and hysteresis damping of rotating shafts. They found 
that with these effects, moments MY and MX in xs and ys planes are 
given by 



( 2 . 2 . 1 ) 

where 7) is the coefficient of internal viscous damping which is 
frequency dependent and tj is the loss factor, representing the 
frequency independent hysieretic damping. 

Applying Newton's law to free— body diagram of a differential 
element shown in Fig. (2.1), one gets 
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Neglecting higher order terms and rearranging, the above 
equations become 

Q' = pAU + cU 

x 

tr + q = pie - i o © 

y x y pox 

Q' = pAV + cV 

y 


M' - Q 
x y 


pi© + I O © 
x pay 


(2.2.3) 


From mechanics of solids, it is known 


Q = k'GA (U' - 6 > 
x y 


Q = k'GA (V' + © ) 

y x 


(2.2.4) 


MX and MY are defined by eqns. (2.2.1). 

Combining eqns. (2.2.1), (2.2.3) and (2.2.4) and arranging 

one gets 


pAU - k'GA (U'~© )' + cU = 0 

y 


pi © -I Cl e — El 
y p o x 


©'' + El lO Y).,+- 

y l ° V | 


=) 

r>J 


-El 77 . , ©" -k'GA (U'-© ) 

V y y 


(2.2.5) 


pAV -k'GA (V'+S ) ' + cV = 0 

x 

.. . r f F/j_j 

pi © +1 Cl e - EI — ©" - EI 0 77 , ,+ 

x p o y U ~J x I o V J oJ 


- EI 7),, ©" +k'GA (V'+© ) = 0 

'V x x 


( 2 . 2 . 6 ) 





15 


Equations (2.2.5) and (2.2.6) are governing differential equations 
in xs and ys planes respectively. These are more general in 
nature and hence are considered from now onwards for analysis. 

2.3 SHAFT WITH UNBALANCE AND NO INTERNAL DAMPING 

Governing differential equations far transverse vibrations of 
rotating shaft, incorporating the effects of shear deformation, 
gyroscopic moment, rotary inertia and unbalance forces are derived 
from Hamilton's principle. 

Consider a differential element of length ds as shown in Fig. 
(2.1). The geometric centre and the centre of gravity of the 
element do not coincide due to presence of unbalance as shown in 
Fig. (2.2). They are related by the equations, 

U = U + e cos & 

9 

V = V + e sin e (2.3.1) 

9 

The kinetic energy of the shaft becomes 

Te = / + V 2 ) ds + J | pi * 6?) ds 

o'" y y o 

£ 1 £ 

+ f ^ I ds + f I deeds (2.3.2) 

J 2 p o J poxy 

o o 

with notations having usual meaning as discussed in section 2.1. 
Plugging eqn. (2.3.1) in (2.3.2) one gets, 

Te = | J jpA (U 2 + V 2 ) + pi (e 2 + e 2 ) 

^ Q v y 

+ 1 £? + 21 £1 6 6 - 2pA U e sin S 6 

p o poxy 
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+ 2pA U e cos 8 8 + e' 




ds 


(2,3.3) 


The strain energy U and variational work done <5W remains as 

e 

discussed in section (2.1) and is given by 


U = f ^ EI (8 2 + e' 2 ) ds 
e J 2 x y 

o J 


/ | k'GA [<u' - e y 

n *- y 


+ (V' + & 


*> 2 ] 


ds 


(2.3.4) 


L L 

6U = J c U «5U ds J c V 6V ds (2.3.5) 

o o 

By Hamilton's principle, 

*2 

6 f (Te - U -W) dt = 0 

and noting that 6 = £1^ , the spin speed of the shaft which is a 

*» 

constant such that 8=0 one gets in the usual manner 

pAU - k'GA (U'-e )' + clt + pAe O 2 cos (Cl t) = 0 

y o o 

pi 8-1 Cl 8 — (EI6' ) ' — k ' GA (U'-e ) = 0 

y p o x y y 

pAV - k'GA (O' +8 )' + cV + pAe CT sin (Cl t) = 0 

X o Q 


pi e +1 Cl 8 -(EI6' )'+k'GA (V'+S ) = o 

x p o y y x 


(2.3.6) 


EI <9' 66 j 1 = 0 , k'GA (U'-e )<5U| l = 0 

x x 'a y ‘o 


EI 8' 68 | L = 0 , k'GA(V'+e )6Vj L = 0 

y y! 0 ’ x >o 


(2.3.7) 
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Equations (2.3.6) gives governing differential equations and 
(2.3.7) gives the corresponding boundary conditions. 

It may be noted that 6 operation on (— pAUe sin 6 6 + pAVe 

cos 6 6) gives terms 




(cos 66U + sin 66V) + pAe (U sin 6 — V cos 6) 66 


Terms with 66 have been ignored in eqns. (2.3.6) as torsional 
vibrations have been neglected in this study. 


2. A NON-DI MENSI ONALI SATI ON 

Let , 


u = 

U 

L 5 

V 

v = r ; 


Cl = 

V*’ 

T* = (p A L 
a o 

4 1/2 

/E I > 

Q Q 

RN 

* p/p o’’ 

YN = E/E ; 

o 

SN = 6/S 

i 

AN 

- S/A i 

IN = I/I 

n 



T = 


(2.4.1) 

Here, L is the overall length of the shaft, p, E, 6, A and I are 

density, Young's modulus, shear modulus, area and moment of 

inertia of the cross-section under consideration, p , E , G , A , 

o o o o 

I are properties at s = 0. Then, 


d n u 

1 

.n 

d u 

d n v 

1 

d v 

ds n 

L n-i 

dS n 5 

ds n 

L n -i 

dS n 

d°U 

L 

d u 

d n V 

L 

,H 

d v 

dt D 

T *n 

dT n ’ 

dt n 

T *n 

dT n 
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d n e 


d n e 

d n e 

. d n e 

y _ 

1 

y . 

x 

1 X 

ds n 

L n 

dS n 

ds 11 

L n dS n 

d n e 


d n e 

d n e 

, d n 

y = 

1 

y 

X 


dt n 

T *n 

dT n 

’ dt" 

T*" df 


Plugging eqns. (2.4.1) and (2.4.2) in eqns. 
and (2.3.6) the non— dimensional ised governing 
differential element becomes, 


(2.4.2) 


(2.2.5), (2.2.6) 

equations far a 


Cl RNM u - SN (u'-Q )' + Cl Dc u = O 

y 


C3 RMN 6 -2 C3 O RMN 6 - Cl RTN 

y * 


1+ *h 




Cl RTN 


^r, H 


Cl RTN r) VN 6" -SN (u'-9 y ) = 0 


Cl RNM v - SN (v'+e )' + Cl Dc u = 0 

y 


C3 RMN e +2 C3 Cl RMN 6 - Cl RTN 
x y 


1+7 Vi 


6 ' 


- Cl RTN 


** 1 +Vj 

1 ,; H 


S" - Cl RTN T ), ... S" +SN iv' +6 ) = 0 

y 'VN x x 


(2.4.3) 


for free vibration analysis and 

Cl RNM u - SN (u'-6 )' + Cl Dc u + RNM e cos (QT) = 0 

y 

C3 RMN 6 -2 C3 Q RMN 6 - Cl RTN S" - SN (u '-6 ) = 0 

y ;< y y 

- 

Cl RNM v - SN (v'+e )' + Cl Dc u + RNM e CT sin (DT) = 0 

X 

C3 RMN e +2 C3 £ RMN 6 - Cl RTN 6" +SN (v'+9 ) = O (2.4.4) 

x y x x 
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far unbalances responses, where , 


RNM = RN * AN; 

RMN = 

RN 

* IN; 

RTN 

= YN 

* IN 

SN = GN * AN; 

Cl = 

Me 

r? 

Re ; 

C3 

= Me 

Re 4 ; 

Me = E /k' S ; 
o a 

Re 2 = 

I 

o 

2 

/A L ; 
o 

^VN 

= 

/T*; 

Dc = cL 4 /E I T 

* 

; « = 

Cl 

<e/L > 





a a 



CHAPTER 3 


TRANSFER FINITE ELEMENT FORMULATIONS 

The governing differential equations for the flexural 
vibrations of the shaft in the non-dimensional ised form are 
derived and given in Chapter II by eqns (2.4.3) for free 
vibration analysis and (2.4.4) for unbalance responses. 


3.1 FINITE ELEMENT FORMULATION FOR FREE VIBRATION ANALYSIS 

The domain of the shaft length is discretised into number of 

finite elements as shown in Fig. (3.1). The typical finite element 

is shown in Fig. (3.2). The variation of u, v, 6 and 6 over a 
* x y 

typical element is approximated as 

(e) t ..u, r . (ne) 

u = [N j -Cu> 

(e) , KI v . r (ne) 

v = [N J -Cv> 

& y = [N y J C6 y > <ne> 


e = [n x j w > <ne> 


6 


(3.1.1) 


6 


where [N^J , [N^ ■> [N , [N X J are shape functions and Iu> 

£v> - {.6 > and IB > are nodal parameters, 

y x 


(ne ) 



i 


(e) 

J (m) 


B) (2) (5) ; C4) < 

2 3 4 5 

H 

S-o 


frl) 

ft (fHl) 

£1 


Fig . 3.1 


Finite Element Discretisation of Shaft. 


iu t 

^ — - 

H 

> iU 2 

ye* 

(e) 

^r©y 

A 

SjJ 

> 6xi 


- ©1:2 

14 h # 


Vt j 

cfH=>0 


d«=H 


Fig- 




Typical Finite Element 
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Substituting eqn. (3.1.1) into eqn. (2.4.3) the residues are 
obtained as 


R U = Cl RNM f N . U u (e> dH - SN f N. U (u' (e) - © <e> ) dH 

J o 1 J 0 1 y 

+ Cl Dc f N. U u <6> dH 

J o 1 


H 


R V = Cl RNM f N. V v (e> dH - SN f" N. V (v' (e) 4- 6 <e VdH 
J x J , i x 


H 

f 

0 


+ Cl Dc f N. V v <6) dH 
O 1 


& H © H & 

R y = C3 RMN f N. y © <e> dH - 2 C3 Q RMN f N. y © <e> dH 

J o 1 y 0 1 x 


y 

H f l+r; H . ©, 

RTN f I 1 N . 

1 


- ci 


v “(e) 
y 6 dH 


+ Cl RTN 


H , -x 6 . 

Jo ( 0 "VN * j== ) N i 


o 

H © 


y ©" (e) dH 
x 


- Cl RTN J r? VN N. 


y g"<e) 


H © 


dH - SN f N- y (u* (e) -© (e> ) dH 

J o 1 y 


x 


H © 


C3 RMN f N. X © <e) dH + 2 C3 (1 RMN f N. 

J o 1 x o 1 


H © 


x q (e) 


dH 


Cl RTN J [ 1 N '* © 

0 477 2 J 1 X 

1+7) H 


(e) 


dH 


H 


Cl RTN J [ r) r .. + 

0 


ri H 6 

2— ] N. 

4 2 J 1 

1+ *H 


VN 

.J ^ 

1+t 'h 

H © .. . . H © 


* ©"< e> dH 

y 


Cl RTN J i) VN N t X ©^ <e> dH + SN J (v 


x - - (e) 


+© ) dH 

x 


(3.1.2) 
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Minimising these residues (3-1.2) by Galerkin method, one 

gets after integration by parts 

H H H 

Cl RNM P N. U ii (e) dH + Cl Dc f N. U u (e) dH + SN f N' U u' (e> dH 

0 1 0 1 J o 1 


H 


SN f N . & 

J o 1 


• u ^ (e) <e> « (e)^ kI u 


dH = SN (lT 


6 N 

y i 


H 


Cl RNM T N. V v 
0 1 


H 


H 


v :-(e) „ r ..v /.(e) r k! ,v , (e> 


dH + Cl Dc f N. 

0 1 


dH + SN f N' v' dH 

J o 1 


H 


+ SN f N'. V e v <e) dH = SN (V 
O 1 X 


' <e) + e (e> > n v 


H 


H e H 6 

C3 RMN f N. y 6 dH - 2 C3 fl RMN f N. y 6 dH 

J o 1 y J o 1 X 


y 

i+n 


h e 


(e ) 


dH 


+ Cl RTN f — If N'. y 6* 

l 2 J 0 1 y 

1 + *H 

r T) w H 6 

- Cl RTN Ci Tj + T N' y 6' dH 

L VN .1 O j J A 1 * 


4 




H 


H 6 H 6 

+ Cl RTN f 7). , h , N' y &' dH - SN f N. y u' dH 
J 0 'VN i y J 0 1 


H 6 


+ SN f N. y e (e) dH 
i y 


o 


?h 


Q ^VN + 


H 


477 2 

l+r?. 


{d RTN [[^=] ^ 
* 1 + r > H 

1 6 ‘ Ce) + r/ uw e' <e> 1 1 n . y 

J ;< VN y J j i 


e ,h 
o 


'H 
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H 6 

r 
0 


H 0 

C3 RMN f N. X 6 <e> dH + 2 C3 Q RMN f N. X 6 (e> dH 
J ,-. i k i y 


+ Cl RTN 


r 1+r; w > H e 

I — ] J N' 

1 477 2 J o 1 


H 

+ Cl RTN f ft n, lkI + — ] f Nt~ X 0' * 

I 'VN j ^ J J 0 1 y 


T)l 


X ©' <e) dH 

X 


H <9 


<e) 


4 


i+'o, 


H 


H & 


h a . . 

+ Cl RTN 7} f N* X Q‘ 6 dH + SN f N. : 

VN J i x J i 

0 O 


H e 

+ SN f N. X Q (e * dH 

J o 1 x 


dH 


, (e ) 

v dH 


e ) 


■ t CI RTN [(p=) 

i+ v 

♦ }}<* 


e ,h 
o 


’H 


(3.1.3) 


Study of the above shows that compatibility of u, v, S^and 

at nodes and completeness in u, v, 6 , 6 , u' , v' , 6' and &' are 

y’ x ’ y x 

required. A linear polynomial i.e. C q elements satisfies these 
conditions. But these lead to ill-conditioning in these kind of 
problems C14,21,243. Hence a cubic polynomial i.e. Cl elements are 


chosen for u 


(e) 


(e) 


a ,e > and e <e> 


a 


The super scripts are 
6 f 

dropped now onwards as [N°J = [N^ = [N ^ J = [N X J - Rearranging 

eqns. <3. 1.3) in matrix form, one gets, 

L N J dH •CuJ <ne) + Cl Dc f CN>L N J dH Cu> <ne) 


Cl RNH J CN> 


+ SN J CN'>1- N 'J dH Cu> (ne> - SN J £N'>L N -1 d H Id j (ne) 
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= SN (u' (e> - 6 <e> -£N> 


H 


H 


Cl RNM J CN> 
0 


L NJ 


dH Cv> 


(ne ) 


H 


+ Cl Dc J CN} L N J dH <v> 


(ne ) 


O 


+ BN J CN' > L N 'J dH £v> <ne> + SN J {N')L N J dH 


(ne ) 


{ P ) < p > 

= SN (v' e +6 CN> 


H 


C3 RMN f -CN>L N J dH id > (ne) - 2 C3 Cl RMN f ^N>i- N J dH CB > (ne ^ 

j a y J a * 


(ne ) 


+ Cl RTN f ] f -CN'} j dH {6 } 

L J 2 J O y 

1+ *h 

- Cl RTN f O T) + 1 f CN' > LN'J d H CBJ 

l VN 4 — 2 ■* 0 


(ne ) 


H 


H H 

+ Cl RTN Th JK , f CN' 3- L N J dH ■£& j (ne - _ r J L J dH id}- 
VN J 0 y J 0 


(m 


+ sn f -CNi L N J dH ce y 

•i . y 


(ne ) 


y 

i+n 


{ ci R ™ [( j==] < e> - ( 


° r; VN + 


'H 


+ r/ VN S y 


B 


H 

£N> 


1 +T), 


=] 


(e) 


t H 


H 
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H 


C3 RMN J £N> 
O 


Ln J 


(ne ) 


H 


dH CB + 2 C3 £5 RMN f im 

K J o 


L N J 


dH W > 

y 


(ne ) 


+ Cl RTN 


f — — — 1 r <n* > L n 'J dH <e > 

^7 0 

( ° * “= ) J^n’>L n ’J dH «y (nei 


H 


(ne ) 


+ Cl RTN 


l+n 


H 


H 


+ Cl RTN r? VN J CN'}^ n 'J dH ie x > <ne> + SN J £N>L N 'J dH £v> ( 


me , 


H , N » 

+ SN f *CN> L -i dH W } 
J X 


(ne > 


= { clR ™ 

^ 1 4-v> 


n . 


e) 


i+T), 


* ^ V ]} 


H 


1+?), 


i H 


£N> 


H 

(3.1.4) 


The above equations are rearranged as 


cm‘ <q> <ne! - [ ccj' - cei* ] tq > < " e > - lki e «,j (ne) = <F> <r ' el 


<3.1.5) 


where 


CM3 


Cl RNM CA1D 
COD 
COD 
COD 


COD COD 

Cl RNM CA1D COD 

COD C3 RMW CA1D 

COD COD 


COD 

COD 

COD 

C3 RMN C1D 
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EG3® = 

E03 

E03 

E03 


E03 



E03 

E03 

EO 3 


E03 



E03 

E03 

E03 

-2 C3 

Q RMN EA1 3 



_ E0:i 

E03 

2 C3 O RMN 

EA1 3 

E03 


EC3 e = 

Cl Dc 

EA13 

E03 

E03 

E03 



E03 

Cl Dc EA1 3 

E03 

E03 



E03 

E03 

Cl RTN 

r) VN CA23 E03 



E03 

E03 

E03 

Cl RTN 

W A2:i - 


SN EA23 
EG3 

-SN EA43 


E03 


E03 

SN EA23 

E03 


-SN LA33 

E03 
A- 


E03 

Sn l A3 j 


Cl RTN 


t^=) 

1+ *h 

SN EA1 3 


EA23 


-Cl RTN 


(" W J=1 


EA23 


H 


SN EA43 Cl RTN 


(“ W ^=] 

1+ V 


EA23 Cl RTN 


P=) 

1+ *H 


EA23 


+ SN EA1 3 


J 


(ne) = L L u J <ne) L v J <ne> L^ v J (ne> L e x J <ne> J 


L U 1 u i u 2 U 2 V 1 “ V 1 V 2 " V i ®yl 6 yl e y 2 *y 2 ®xl ~®k 1 S x 2 ~ 6 x 2 J 


F> <ne> = {_ _SLf i 0 SU o O -SVj O SV 2 O — MY 1 O MY 2 O -MXj O MX 2 O J 


(3.1.6) 
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where 


CA13 

= 

n 

J 

0 

{N5 L N J 

dh 

CA23 

— 

H 

J 


dH 



0 



CA33 

= 

H 

J 

tN'jLNj 

dH 



0 



CA43 

— 

H 

J 

L N 'J 

dH 



0 



SU 1 

= 

SN 

(u Me) - 

JY 

y 1 1 

su 2 

= 

SN 

, ' (e) 

(u - 

^ e \ 

y*2 

SV 1 

= 

SN 

, ' <e) ^ 

( v + 

e (e, L 

X 1 1 

sv 

= 

SN 

, ' <e> ^ 

( v + 

x *2 


MY i 

MY 2 

MX X 

MX ? 


| Cl RTN 
| Cl RTN 
| Cl RTN 
| Cl RTN 



^ ^VN + 


1 ft* ^ ® ^ J_ yw £i* C e ) 

J x * "vN y 


0 "VN + 


) e i <el+ "vn «; 




° r VN + 


= } e y <6 '- 


S' (e> 


Q ^VN + 


) e ;‘ e>+ K 


* <e) ] 


(3.1.7) 
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3. a FINITE ELEMENT FORMULATIONS FOR UNBALANCE RESPONSES 

Discretising the shaft, approximating the variations of u 


(e) 


(e ) ^ (e ) 

v , & 

y 


and 6 


(e) 


and applying Galerkin method to the 


residues obtained as discussed in the previous section, eqns. 
(2.4-4) that govern the unbalance response become in the matrix 
form as 


CM] {q> 


q> (ne) + [cC3 e + C63 e ]£q> (ne) 


+tK3 e {q> ne) 


{Q3 <ne> +£F> <ne> 


(3.2.1) 


where 


= 

Cl RNM EA13 

CO 3 

E03 

E03 


E03 

Cl RNM EA13 

E03 

E03 


E03 

i 

E03 

C3 RMN EA13 

E03 


j EQ3 

E03 

E03 

C3 RMN 


EG3 


sr 

E03 

E03 

E03 

E03 


E03 

E03 

E03 

EO 3 


E03 

E03 

E U J Cw* 

Q RMN 


E03 

E03 

2 C3 a RMN EA1 3 

E03 


= 

Cl Dc EA13 

E03 

E03 

E03 


C03 

Cl Dc EA13 

E03 

E03 


E03 

E03 

E03 

E03 


E03 

E03 

E03 

E03 



3i 


CKD e = 


SN EA 2 D COD 

COD SN CA 2 D 

-SN EA 4 D COD 

COD SN CA 4 D 


-SN CA 3 D 
COD 

Cl RTN CA 2 D 
+SN CA 1 D 

COD 


COD 

SN CA 3 D 
COD 

Cl RTN CA 2 D 
+SN CA 1 D 


L q ] <ne) = L j <ne) LV J <ne) L a v J <ne> LA j <ne> J 


[Q> 


(ne) 


tFi <ne> 


where 


L U 1 u 2 u 2 V 1 ~ V i V 2 ~ V 1 °yl °yl w y2 w y2 ~xl “xl ~x2 “x2 


H H 

[_RMNj e n 2 cos OT* [Nj RMNj o e fi 2 sin OT* [Nj [Oj [Ojj 


e „ e* . e „ ©*.„ 

T 


j_ _-SU 1 0 SU, 0 -SVj 0 SV 2 0 -MY 1 0 MY 2 0 -MXj. O MX 2 O J 

( 3 . 2 . 2 ) 


T 


H 


CA 1 D = 


CA 2 D = 


CA3D 


CA 4 D = 


sir 


SU, 


SV 


L«J 

fN' > L N J dH 

J CN> L N 'J dH 
o 


J £N> C_ ■' _J dH 
0 

H 

J CN'D 

0 

H 

J £N' D 1 — " " — 1 dH 
0 


sn <u' <e) - e y <e) )| 1 


SN (u' (e) - e y <e> >j 2 
SN ♦ »„<•»>!, 
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SV„ = SN (V* (e) + e (e) ! 

- x <2 

MY. = Cl RTN 6 ’ (e) I 

1 y 1 1 

MY = Cl RTN B* <e> | 

2 y »2 

MX, = Cl RTN (e> ! 

1 x 1 1 

MX = Cl RTN 0' (e) { (3.2.3) 

3.3 FORMULATION OF FINITE ELEMENT TRANSFER MATRICES 

Equations (3.1.5) represents the flexural vibrations of the 
beam for free vibration analysis. Eqn . (3.2.1) are the governing 

equations for unbalance responses of the rotating shaft. Finite 
element transfer matrix formulations of both cases are discussed 
here . 


3.3.1 Free Vibration Analysis 

Assuming the solution to be harmonic, the nodal displacements 
can be written as 


(ne ) r , (ne ) XT 
= -Cp> e 


(3.3.1) 


where p are magnitudes and X is the natural frequencies of whirl 
(both non— dimensionalized and complex) C193, as 
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y y 

Here will give the magnitude of nodal displacement and 
tan * Pj/P^ represents the phase, a is the growth factor and to is 
the natural frequency of whirl. At critical speeds, the rotating 
shaft speed Cl is equal to to for forward and —to for backward whirl 
speeds. 


Introducing eqn. (3.3.1) in eqn. (3.2.5) one gets 


jx 2 CM3 (e) +X. jce3 <e> +CC3 <e) j+CKD <e) jcp> 


(ne) _ ^^(ne) 


(3.3.3) 


Here , 


[F5 <ne> - (F R + iF I > <ne> 


(3.3.4. ) 


as similar to *Cp> in eqn. (3.3.2). 

Substituting eqns . (3.3.2) and (3.3.4) in eqn. (3. 3.« 3) 

one gets 


K- 


-to^) CM3 + <y jjlGD + CCD j + CK3 j + i £ (2 


eto) CM3 + 


id 

CSD + CCD 

|] 

{ Pr + ip i} 

= 

{ f k + iF I } 




\ ■r 




(3.3.5) 


Letting 

(o^-cS 2 } CM3 + a CCS3 + CC33 + CK3 
CZ&iti) CM3 + 6j CCG3 + CC3 3 = CAI3 
one gets 

f CARD + i CAI3 


CARD and 


]{ P R + iP I } { F R + lF I } 


(3.3.6) 


Separating real and imaginary parts and rearranging, one gets 



AR 1,1 

A1 l,l 

- AI 1,1 

AR 1,1 

AR 1 ,2 
AI 1,2 

A I . ^ - - * 

1 

AR 1,2 ••• 

AR i,e 

AI 1,8 

- AI 1,8 
AR 1 ,8 

AR 1,9 

AI 1,9 

- AI 1,9 
AR 1,9 ••• 

AR X 

AI 1. 

,16 

f 16 

- AI i, 

AR 1 

,16 

,16 

fuR, 

ul 

uR 













:ul 













| : uR. 













1 ^2 


^ sim 1 Y 


-sui 


i ! 


o j 

I 

0 ! 

5UFU 

SUI, 


AR 8,1 

I 

3 > 

►H 

CD 

AR 8,2 

- ftI 8,2 --- 

AR 6,8 

1 

T > 

n 

CD 

CD 

AR 6,9 

AI 8,9 

" • * 

AR B ,16 

AI 8,1 

AR 8,1 

AI 8,2 

AR 8,2 

AI 8,8 

AR 8,8 

AI 6,9 

AR B ,9 

* * ’ 

AI 8,16 

AR 9,1 

- AI 9,1 

AR 9,2 

- A 1 9,2 --- 

AR 9 ,B 

" A 1 9,8 

AR 9, 9 

_ AI 9,9 

* * * 

AR 9,16 

! AI 9 , 1 

AR 9,1 

AI 9,2 

AR 9,2 --- 

AI 9,8 

AR 9 ,e 

AI 9,9 

AR 9,9 


AI 9,16 


o 

o 


n 


-MYR, 


-MY I, 


5UR„ 


2 ■! 

| SUI^ 


AR ! 6,1 - AI 16,1 AR 16,2 " AI 16 ,: 


AR 


-A I w 


16,8 16,8 AF 16,9 "‘ 16,9 


— AI , 


< AR 16,1 A 1 16,2 AR 16 ,: 


AI, 


16,8 AR 16,8 H ‘ 1C .9 - 16,9 


A I , 


AR, , 


16, 1 


16 , 16 

AI 16, 16 ! 

- e # R « 

X d 

0 

16,16 

AR 16,16 

- e * 

1 x Y l 

V 

0 


( 3 - 3 * 7 ) 
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Grouping the variables of node 1 and node 2 separately, the 
above equation can be re— written as 


AB 1,1 

AB 1 ,2 

m m m 

AB 1,16 

AB !, 

17 

AB 2,1 

AB 2,2 

mum 

AB 2, 16 

AB 2, 

17 


AB 16,1 

AB 16,2 

9 ** m 

AB 16,16 

AB ld,,17 

m m m 

AB 16,32 


6' I, 
x 1 

-_z_ 

0 

sufC,^” 

AB 17,1 

AB 17,2 


AB 17, 16 

AB 17,17 

« m m 

AB 17,32 

' 1 

! 

| uR 2 ( 




• - - 



9 m m 

j 

1 

i ul 

i 

! 

i 

s 

1 

1 

SUI 2 

1 

mm m m m 

AB 32,1 

m m m m m. m 

AB 32,2 


AB 32,16 

| AB 32,17 

m m m 

AB 32 , 32 

! ^ R 2 

E > 


° ; 


AB 


AB 


1,32 

1,32 


f uRj ^ 

UI 1 

vRl 

vl. 


e r h 

y i 

& i, i 
y i : 

e r . 

X 1 


f-SUR ^ 
-SUI. 


-SVR. 


-SVI 


-MYR. 


—MY I . 


-MXR, 


(3.3.8) 

Equation (3.3.8) is partitioned as shown and is re-written 


as 


CAB13 

18x16 

EAB33 

16x16 


EAB23 

16x16 

EAB43 

16x16 



f 

kp> x 


r 

£-F>i 


16x1 


16x1 

1 

£p> 2 
i 16x1 

r = " 

cf> 2 

16x1 


L J 


L J 


(3.3.9) 
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Further rearranging eqn. (3.3.9) one gets 


{ 


< P> 2 1 


-CAB23 1 CAB13 

- CAB23 1 


r % 

<F> 2 j 

► = 

£AB33-CAB43tAB23 _1 EAB13 

- CAB43CAB23 -1 


1-..1 


or 


CW> 2 = ET3 1 


(3.3.10) 


Here ETl^ is the transfer finite element matrix relating the 
variables at node 2 to those at mode 1. 


3.3.2 Unbalance Response Analysis 

2 

The unbalance present will induce a force given by <eQ so 
that the whirl speed becomes the spin speed of the shaft. Hence 
at steady state, the response will be given by 

= £p> e (3.3.11) 

where p are complex as discussed in the previous section, i.e. 

o = o + i o _ (3.3.12) 

F F R Fj 

Introducing eqns. (3.3.11) and (3.3.12) in eqn. (3.3.1), one 


gets 


^-a 2 EM3 <e> + iO^CG3 (e) + CCD <e> j + CKD te> J{p R +iPj> 

= £Q D +iQ T > <ne> + (F. + iF } (nel 
R I R I 


(ne > 


(3.3. 13> 


Separating real and imaginary parts as discussed in Sec. 
and re-arranging, one gets. 


(3.3.1 ) 
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Shifting the unbalance farce vector into the coefficient 
matrix the above equation are further be rearranged as E53. 

Eqn. (3.3.6) is re-arranged and obtained as : 
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AB 


1,1 


AB 


1 


AB 


2,1 


AB, 


AB 16,1 AB 16,2 


ab 17,i ab 17,: 


AB 


AB 


1,16 

2,16 


AB 1,17 AB 1,32 

AB 2,17 *'* AB 1 ,32 


AB 


16,16 


. AB 


AB 


16,17 


AB 


16,32 


AB 

17,16 | 17,17 


... AB 


17,32 


AB 32,1 AB 32,2 


AB 32,16| 32,17 "* AB 32,32 




- Q 1 

/ 

uRj ’ 

1 0 


Ul l 

-% 


vRi 

0 


Vl l 

0 


Vl 

0 


6 I , 
y i 

0 


6 R 

X 1 




1 0 


e ' i 

2LJ) 

1 _Q 3 


uR 2 

0 


uI 2 | 

-s 


VR 2 

m m m 



1 ° 
a 

. i 

*^2 
l / 


(3.3.15) 


The above equations are partitioned as shown 


and rewritten as 


EAB13 

CAB23 

CAB5> 

1 6x 1 6 

16x16 

16x1 

CAB33 

CAB43 

CAB6> 

16x16 

16x16 

16x1 

L°J 

L°J 

1 

1x16 

1x16 



f ^ 

Cp>i 


— 

iH 

a"S 

U- 

1 

J 

16x1 


16x1 

Cp > 2 

h = 

C f > 2 f 

16x1 


16x1 

1 


1 



L. / 


(3.3.16) 
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£a>: ( 3 .3 6; is re -arranges 


5ined as. 


1-EAB23 i CABi 3 


-CAB23 


-CAB23 -CABS > 


j CAB33— CAB43EAB23 EAB13 -EAB43 EAB23 1 -EAB&}— EAB43 EAB23 1 EAB5> 


(3.3. 17 > 


33x1 


where CT3 ^ is the transfer finite element matrix that relates 
variables at node 2 to those at node 1. 


3.4 APPLYING BOUNDARY CONDITIONS 

Transfer finite element matrices for each element is 
evaluated as discussed and multiplied over entire shaft length as 


IW> = CT3 . CT3 , 
n n-1 n-^ 

CW> = COT 3 -CMl 
n 1 


CT3 7 CT3 j CW> 1 or 


(3.4.1 ) 


where C0T3 is overall transfer finite element matrix relating 

variables at node n to those at node 1. In eqn. (3.4.1) the 

nodal variables at node n and node 1 are given by 

|Wl = luR ul vR vl 6 R 6 1 6 R 6 1 uR' ul' vR' vl' 6 R' 

L J n L u, ' n u n n nynynxnxn n n n nyn 

6 I' 6 R' 6 I' SUR SUI SVR SVI MYR MY I MXR MX I 
ynxnxn n n n n n n n n 

0000000 Oj 
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L W J ! = L uR l Ul l vR l vI i e y R 1 Vl S X R 1 Vl uR l Ul i vR i Vl i e y R i 

© Ij © X I' SURj SUI 1 SVR 1 SVI 1 MYR 1 MY^ MXRj MXI 1 

OOOOOOO Oj (3.4.2) 

Method of applying boundary conditions corresponding to a 
simply supported rigid bearing and elastic hydrodynamics bearings 
is discussed below. Other types of boundary conditions can be 
incorporated in a similar manner. 


3.4.1 Simply Supported Rigid Bearings at Both Ends 

The imaginary part of the complex variable at any node 
indicates the phase shift with respect to vriables at node 1. 
Hence the imaginary part at node 1 are all individually zeros. 
Hence from eqn . (3.4.2) -CW1 ^ reduces to 

[WJj = lu*! e y Ri ®X R 1 uR i uR i e y R l S x R i SUR 1 SVR i MYR 1 MXR 1 J 

(3.4.3) 

For simply supported rigid bearings, the boundary conditions 

are 

u = (uR + iul ) = v = (vR + ivl > = MY = (MYR +iMYI > 
n nnn nn n nn 

= MX = (MXR +iMXI > = uR. = vR = MYR = MXR = 0 (3.4.4) 

n n n 1 1 1 1 

Combining eqns . (3.4.1), (3.4.2), (3.4.3) and (3.4.4) one 

gets 

tw> = COT 3 -CWK (3.4.5) 

n 1 

16x1 16X8 8x1 
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Here -CW> ^ corresponds to vanishing variables at nodoe n and *CW> ^ 
corresponds to nonvanishing variables at node 1, which are 

[w I = L uR r, uI „ vR vI MYR MY I MXR MX I OOOOOOOOl 
L J n *- n n n n n n n n J 

L W Jl * l 6 y R ! e x R ! uR i vR i e y R j V4 SUR X SVRj (3.4.6) 

Equation (3.4.6a> corresponds to a combination of real and 
imaginary parts of variables at node n which can be seperated 
into two matrices as 

iWR 1 = luR vR MYR MXR 0 0 0 0 I 

LJ n L nnnn -* 

|WI 1 = |ul vl MY I MX I 0 0 O 0| (3.4.7) 

LJ n L nn n n J 

Combining eqns. (3.4.5), (3.4.6), and (3.4.7) one gets 

CWR> = COT _.3 -£W> 
n R 1 

8 xl 8x8 8x1 

CWI> = COT_ 3 -CW> (3.4.8) 

n II 

8 xl 8x8 8x1 


3. 4. 2 Elastic Hydrodynamic Bearing Support at Both Ends 

When the end supports are elastic bearings, then the boundary 
conditions in non-diamensional ised form become C 15, 173 
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where the bearing stiffnesses and dampings are 
diamensional ised as 


k 


XX 


k L 

Cl and C 

El xx 

o o 


Cl 


c L 

XX 


E I T 
o o 


Assuming harmonic solution of u, v, 0 and 0 as 

x y 

XT XT ^ ^ XT ^ ^ XT 

u=ue , v = ve , 0 = 0 e ,0 = 8 e , 

y y x x 

where X = cy + ico and re-arranging eqns. (3.4.9) into real 
imaginary parts, one gets 


SUF^ = SN 


f(K , 

[ xx 1 


oC , ) uR, - a) C , ul 
xxl 1 xxl 1 


+ (K , + oC , > vR, - ci C , 
xyl xyl 1 xyl 


-i] 

SUI 1 = SN [ “ c xxl uR l + (K xxl + OC xxi ) ul l 

'xyl' ’*1 ] 




+ w C , 
xyl 

vR x + 

(K + oC 

xyl : 

SVR t 

= SN 

K*i * 

aC 4 ) 
yxl 

uR^ - ci C 



+ (K • 

yyi 

*■ oC , 

yyi 

3 

1 

•jH 

cc 

> 

SVIj^ 

= SN 

[ “ c y*i 

uRj + 

(K , + OC 
yx 1 



+ “ c yyi 

vP^ + 

(K , + oC. 

yyi 


yyi 

yxl* 

yyi } 


■3 


] 


SUR n * SN ♦ oC xxn ) uR n + * C x)<n ul n 

- <K *yn * ^xyn 1 uR n * “ C xyn vI n ] 


non— 


and 
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SUI = SN 
n 


SVR = 
n 


SVI 


MYR. = 


[-CD C 

uR - 

(K 
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oC ) 

ul 

[_ xxn 

n 

xxn 


xxn 

n 

- CD C 

vR - 

(K 

4- 

oC ) 

vI 1 

xyn 

n 

xyn 


xyn 

n J 

[-(K 

+ cC 

) uR 

+ 

03 C 

ul 

L yxn 

yxn 

n 


yxn 

n 

(K 

+ oC 

) vR 

4- 

63 C 

vl 1 

yyn 

yyn 

n 


yyn 

n J 

I CD C 

uR - 

(K 

4- 

oC ) 

ul 

L yxn 

n 

yxn 


yxn 

n 

" CD C 

vR - 

(K 

4- 

oC ) 

VI 1 

yyn 

n 

yyn 


yyn 

n J 

= MXRj = MX I 

= MYR 

n 

= MY I 

n 

= MXR 


MX I = O 
n 


(3.4.10) 


At bath ends the nodal farce variables are either zeros or a 
function of nodal displacements and bearing properties. Hence 
these can be transferred to corresponding coefficients in over all 
transfer matrix. It is to be noted that imaginary parts of 
variables at node 1 are zeros. Therefore, combining eqns. 
(3.4.1), (3.4.2) and (3.4.10) one gets 

-CW> = COT 3 -CW3 _ (3.4.11) 

n 1 

16x1 16 xS 8 x 1 

Here £W> corresponds to vanishing variables at node n and tW> 
n n 

corresponds to non— vanishing variables at node 1, which are 


L»J„ - 

Wt - 


LSUR, SUI p SVR n SVI n MYR n MYI„ MXR r MXI n 

L uR 1 YR t e y R, e x R i vr; e y R ; j 


ooooooo oj 

(3.4.12) 
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Equation <3.4. 12a) can be separated into two matrices 
corresponding to real and imaginary parts as 


L«"J n - 

ISUR SVR MYR MXR 

L n n n n 

0 0 0 Oj 


L WI Jn " 

ISUI SVI MY I MX I 

L n n n n 

0 0 0 Oj 

<3.4.13) 


Combining eqns. <3.4.11), <3.4.12) and <3.4.13) one gets 

CWR> = COT 3 

n R 1 

8 xl 8x8 8x1 

CWI> = COT 3 CW>, <3.4.14) 

nil 

8 xl 8x8 8x1 

3.5 APPLYING BOUNDARY CONDITIONS FOR UNBALANCE RESPONSES 

As discussed in Sec. <3. 4), one gets over all transfer finite 
element matrix as 

CW> 7 = CTlj CW> 1 <3.5.1) 

33x1 33x33 33x1 

This can be separated and rewritten as 


f ^ 

CWJ 

n 

32x1 


CAB13 

32x32 

CEF> 

32x1 


r 1 

32x1 

’1 

1 


L°J 

1 

1 

r 

1 

u J 

! 

1 

1x32 



L - 


where CW> , CW> „ are nodal variables as discussed in eqns. <3.4.2) 
n T 1 

CEF> corresponding to the external force matrix, due to mass 
unbalance . 

Boundary conditions are now incorporated as discussed in Sec. 


<3.4.1) and <3.4.2) to obtain 
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(3.5.3) 


Here again -CWJ^ corresponds to vanishing variables at node n 

and £W> to non— vanishing variables at node 1. As discussed in 
n 

Bee. (3.5), eqn. (3.5.3) can be separated corresponding to real 
and imaginary terms as 


xWR> 

n 

- co v 

{W} + -CEF > 

1 K 

= -C0> 


8 xl 

8x8 

8x1 

8x1 



CWI> 

n 

= COTjl 

(W) + CEFjl 

= C0> 

(3.5.4) 

8 xl 

8x8 

8x1 

8x1 





CHAPTER IV 


RESULTS AND DISCUSSIONS 


Computer program in FORTRAN 77 have been developed to obtain 
natural frequencies, critical speeds and unbalance responses. 
Double precisian was needed to obtain these results. 

4.1 SOLUTION TECHNIQUES 

4.1.1 Free Vibration Analysis 

Assume a starting value of ct and w (say o^and ) and obtain 

eqns. (3.4.14). If assumed value of a and o (a , oa ) are the 

o o 

correct roots then 
D1 = |0T_ 


D2 = jOT 


= 0 
= o 


(4.1.1) 


Since a and w are assumed values, eqns. (4.1.1) are not 
o o 

likely to be satisfied. Correct values of ct and w can be obtained 
by using Newton— Raphson' s iteration method 119,223. Thus 


r *> 



► = 




L 


8D1 


0D1 

did 


-,-1 


da 

dD2 &D2 


L da 


d<d J 


-D1 ict t id ) 
o o 

-D2 ia ,<d) 
a o 


(4.1.2) 


so that 



The iterative process is repeated till da and du falls below 


a pre-defined small value. Derivatives in eqn . (4.1.2) are 

evaluated as 


dD l 

D1 

ict ) - 

o ’ o 

D1 

(a ,63 ) 

o’ o 

da 


6a 




D1 

(a ,63 +<563 > - 

o o 

D1 

ia ,63 ) 

o o 

<?63 


{£’63 



<?D2 

D2 

(a +£a, 63 ) - 
o o 

D2 

( a ,63 ) 

o o 

da 


6a 



£D2 

D2 

(a +6t*) - 

cr o 

D2 

{a ,63 ) 

o o 


<563 


(4.1 .4) 


Thus to evaluate these derivatives, evaluations of D1 and D2 

at (<? , u ), (a +&a. 63 ) and ia ,63 +6oj) are required, 

o’o o o oo 

In Newton-Ralphson method, convergence of the results is 

dependent on a good choice of starting value of c and 63^. This 

is a serious problem when damping either internal or external is 

present. Experimentation in this work has shown that this is best 

done, first by evaluating the undamped natural frequency. For this 

a was assumed to be of the order of 10 6 and 63 chosen 
o 


arbitrarily. Now, 1)4 of the damping is introduced. 


For this. 


starting values of a and 63^ are chosen at 2 to 3 and 90)4 of 


undamped natural frequency respectively. The roots O' and 


63 thus 
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obtained serve as the starting point for the next step of damping 
which was chosen at 10% whose roots were chosen to be starting 
value for 100% damping. The above is summarised in Table 4.1. 

It was also found that £& and at 0.1% of the present value 

°o anc * 9* ves fast converging results for this kind of 
problems. 

4. 1 . 2 Response Analysis 

Eqns. (3.5.4) are re-written as 

_ 0T r] Wi + { EF r} ' { WR }n - {«} 

= or 

Hi - - Hr] { ef r} ' -[° T i] 1 HO (4 ' 1 - 5 ’ 

Either of the two of the above equations can be evaluated to 
obtain non— vanishing variables at node 1. Both were found to give 
identical results as expected. 

Thus, having found variables at node 1, variables at other 
nodes can be evaluated using eqn. (3.3.8) as 

CW> = CT3 , -CW> (4.1.6) 

n n-1 n-1 

4.2 SHAFTS WITHOUT DAMPING 

Critical speeds in forward precession were obtained solving 
eqn. (3.4.1) with all damping values taken zero for a uniform 
steel shaft whose material constant (Me) is 2.9. Results were 
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Table 4.1: 


Step 


1 


Procedural Steps 


Damp ing 


0 


for Evaluating Complex Frequencies. 


Starting Value Converged Value 


Cf LA Ct LA 


10 ** arbitrary 0 la^ 


1 % 


-2 to —3 0.9 la 

i 


a 


2 


LA 


2 


3 ioy* 


ct 

2 


<A 


2 


ct 

3 



4 iooy 


a la ct (A 

3 3 4 4 

( roots ) 
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obtained far rotary inertia coefficient (Re) equal to 0.05 and 
0.01^5 which represents a shaft having length to diameter ratio 5 
and 20 respectively for h inged— h inged and cantilever boundary 
conditions. Also results were obtained for different end 
conditions in two planes. 4 , 6 and 8 elements were used and 
results are compared with £143. These are given in Tables 2, 3 
and 4; these match very well with the earlier results £143. 

Reference £143 found that using 8 elements i.e. 36 degrees of 
freedom gave accurate results when two second order equations were 
used. In the present work using four second order equations, it 
was found 4 elements (amounting to 40 degrees of freedom) are 
sufficient to give accurate results with much reduced memory 
requirement. Hence now onwards, four elements are considered for 
study in this work. 

4.3 SHAFTS CARRYING DISK WITHOUT DAMPING 

Critical speeds are obtained for a steel shaft of 1000 mm 
length and 40 mm diameter with hinged-hinged supports carrying a 
disk of 30 mm thickness and 90 mm dia placed at 250 mm from the 
left end, (Fig. 4.1). The disk is also considered flexible. 
Effect of shaft mass to disk mass is also studied. Range of shaft 
mass is obtained by successively reducing the density of the shaft 
maintaining all other parameters constant. The results are given 
in Table 5. They match very well with earlier results £143. 

It is seen that there is an error of 158% for a shaft to disk 
mass ratio of 6.584 and 30% for shaft to disk mass ratio of 0.844. 



Study of Critical Speeds for Hinged— Hinged . 
Shaft . 




Re = 0.05 

Me = 2.9 


Mode 


No. of Elements 


Mahendrakar 

r •§ ii t 

4 

6 

8 

L 14 J 

8 Elements 

i 

9.6384 

9.6384 

9.6383 

9.6383 

2 

35 . 8775 

35.8775 

35.8774 

35.8774 

T 

vJ 

72.1430 

72. 1430 

72.1429 

72. 1429 

4 

112.2963 

112.2960 

112.2959 

112.2958 

5 

153.2134 

153.2124 

153.2123 

153.2122 

Table 

4*3: Study of Critical Speeds for 

Cantilever Shaft* 



Re = 0.0125 

Me = 2.9 


Mode 


No. of Elements 


Mahendrakar 

C 143 

8 Elements 

4 

6 

8 


3.5134 

3.5128 

3.5123 

21.9454 

21 .9404 

21.9363 

61.0262 

61.0102 

61.0101 

118.4774 

118.4604 

118.4583 

193.5804 

193.5723 

193.5702 


3.5123 

21.9363 

61.01 

118.46 


5 


193.5702 
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Table 4.4: Study of Critical Speeds with different support 
conditions in two planes 

Re = 0.05 Me = 2.9 


Hinged 


Hinged 



Fixed 


Fixed 


Mode 

No. 

of Elements 


Mahendrakar 

L 143 

8 Elements 

4 

6 

8 

i 

9. 6408 

9.6402 

9.6383 

9.6383 

2 

19.3448 

19.3423 

19.3413 

19.3401 

3 

38.8962 

35.8836 

35.8800 

35.8800 

4 

47.0583 

47.0453 

47.0403 

47.0400 

5 

72.1563 

72.1402 

72. 1483 

72.1400 

6 

81 .3672 

81.3500 

81.3482 

81.3500 
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Thus , one should not neglect the shat t mass, as error can be 
considerable even in first mode when shaft and disk mass are 
comparable . 

4. 4 SHAFTS WITH DAMPING AND ELASTIC BEARING SUPPORTS 

To study effects of internal viscous and hysteresis damping, 

eqns. (3.4.1) were solved. A steel shaft with Re = 0.02 (i.e. 

length to diameter ratio 12.5) studied by [12, 141 with bearing 

supports having uncoupled stiffness and damping characteristics 

(i.e. K =K = C =C =0) was chosen. Each case is listed 
xy yx yx xy 

below . 

Case 1: Viscous damping, isotropic undamped bearings 

Bearings are undamped and have isotropic bearing conditions 

of K = K = 33.14 ( = 17.5 MN/m) 
xx yy 

Viscous internal damping is considered with r, VN = 0.01621 

(0.0002s) . Natural frequencies and growth factor are given in 
Table 6, as a function of shaft spin speed fl. The first critical 
speed occurs at 6.219 (4805 rpm) and second at 13.0043 (10048 rpm) 
in forward made of precision. The critical speeds observed by 
C12, 143 were at 4980 and 10470 respectively, which are higher 
than the present value as C12, 143 had neglected effect of shear 
deformation . 

The viscous internal damping begins the destabilising effect 
only after a certain spin speed is reached. It was found that the 
first mode is stable up to the first forward critical speed and 
instability sets in then. The second mode becomes unstable at the 
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Table 4,6: Natural Frequencies of Shafts with Internal Viscous 
Damping and Isotropic Undamped Bearings 


K 

XX 

K = 33.14 C = C 

yy XX yy 

0.0 

= 0.01621 

Spin Speed 
rpm 

Growth 

Factor 

Natural 

Frequency 

Present 

Nelson E123 

fNon-dimensionaiised' 

Present 

Nelson C 123 

(From graph srd 
Nor'- Sur er s ’si 

CaD Forward Precession 





-.1124 

-.1321 

6.232 

6.4 

20 OQ 

-.1102 

-.1116 

13.462 

13.6 


-.0643 

_ 

6.221 

6.4 

3000 

-.0765 

— 

13.323 

13.6 


.0024 

.0010 

6.201 

6.4 

4950 

-.0327 

- . 0605 

13.131 

13.6 


.1426 

.1303 

6.172 

6.4 

10000 

-.0023 

-.0053 

13.001 

13.6 


.2026 

.1909 

6.163 

6.4 

12000 

. 0420 

.0220 

12.963 

13.6 

CbD Backward 

Precession 




-.1288 

-.1322 

6.424 

6.4 

2000 

-.1114 

-.1116 

13.673 

13.6 


-.2122 



6.462 

6.4 

3000 

-.2024 


13.683 

13.6 


-.3137 


6.563 

6.4 

4950 

-.2983 

-.1552 

13.730 

13.6 


-.3347 

-.3830 

6.731 

6.4 

10000 

-.1785 

-.1925 

13.824 

13.6 


-.3243 

-.4394 

6.820 

6.4 

12000 

-.1988 

-.2074 

13.900 

13.6 
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second critical speed. This is in agreement with observations 
made by C12, 143. Both these modes remained unstable from their 
critical speeds upto the speed range studied. 

Common observations for all the cases are listed separately 
after Case 5. 

Case 2: Viscous damping, isotropic damped bearings 

The system studied now is Case 1 with isotropic bearing 
dampings, C^ = C = 0.269 (17.5 KN.s/m). The results are shown 
in Table 7 . The first mode is now stable up to a spin speed of 
11.23 (11.6 by ref. C 12 , 143 ) and second mode up to the entire 
speed range studied. This matches with observations made by 
C 12, 143 . 

Case 3: Viscous damping, anisotropic undamped bearings 

The shaft has a viscous internal damping = 0.01621 and is 

supported in two similar undamped anisotropic bearings whose 

stiffness are K = 33.14 and K = 66.32. 

x:< yy 

The results are shown in Table 8. The system is stable for 
the entire speed range studied which is in agreement with 
observations of C 12 , 143. 

Case 4: Hysteresis damping, isotropic undamped bearings 

The effect of hysteresis internal damping is studied in this 
case. The system is Case 1, with = 0 and loss factor = 

0.0002. Results are given in Table 9. It is seen that both 

forward precession modes are unstable for the entire speed range 
studied. Same was observed by C12, 143. 
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Table 4.7: Natural Frequencies of Shafts with Internal Viscous 
Damping and Damped Istropic Bearings 


K 


xx 


= K 


yy 


33. 14 


x:< 


= C = 0.269 

yy 


^VN = °- 0i621 


Spin Speed 

Growth 

Factor 

Natural 

Frequency 

rpm 

Present 

Nelson C123 

(KJor> -di’merisional is ed] 

Present 

Nelson 1121 

( From arBph and 

Mon -dim &i S /'ons \i sedj 

C sO Forward Precession 

-.2111 

-.2329 

6.293 

6.35 

zwu 

-.7903 

-.7971 

13.11 

13.55 


-.1027 


6.20 

6.35 

3000 

-.7329 

- 

13.005 

13.55 

4950 

-.0943 

-.1002 

6.185 

6.35 

-.7432 

-.7662 

12.897 

13.55 

8000 

-.0322 

-.0552 

6.131 

6.3 

- . 6838 

-.7133 

12.702 

13.5 

10000 

.0428 

.0286 

6.108 

6.2 

-.6117 

-.6746 

12.543 

13.45 

Cb) Backward 

Precession 

-.2124 -.2329 

6.36 

6.35 

2000 

-.6913 

-.7971 

13.602 

13.55 


- . 2326 

_ 

6.355 

6.35 

3000 

-.7124 

- 

13.676 

13.55 


-.3431 



6.464 

6.37 

4950 

- . 7324 

-.7662 

13.682 

13.55 


-.3932 

- . 4260 

6.625 

6.4 

8000 

-.7634 

-.8680 

13.743 

13.6 


- . 4343 

-.4853 

6.700 

6.5 

10000 

- . 7843 

-.8873 

13.794 

13.65 
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Table 4.8: Natural Frequencies of Shafts with Internal Viscous 
Damping and Undamped Anistropic Bearings 


K 


xx 


14 


yy 


= 66.32 C 


xx 


= C = 0.0 n 

yy 


VN 


= 0.01621 


Spin Speed 

Growth 

Factor 

Natural 

Frequency 

irpfn — 

Present 

Nelson C123 

Present 

Nelson C12 1 



(Ki'on-dirnensionsiised} 


(From graph and 
Non -dimensiora J/s en 

C a.) Forward 

Precession 




500 

-.1241 

-.1318 

7.423 

7.6 

-.1324 

-.1542 

17.921 

18.5 

2500 

-.1164 

-.1297 

7.420 

7.6 

-.1250 

-.1518 

17.902 

18.5 

5000 

-. 1040 

-.1226 

7.416 

7.6 

-.1194 

-.1434 

17.882 

18.5 

8000 

- . 0943 

-.1114 

7.384 

7.6 

-.1046 

-. 1303 

17.692 

18.5 

10000 

- . 0648 

-.0960 

7.301 

7.6 

- . 0947 

-.1123 

17.502 

18.5 

C fcD Backward 

Precession 



500 

- . 2583 

-.2725 

6.452 

6.3 

-.2128 

-.2329 

13.851 

13.5 

2500 

-.2686 

- . 2737 

6.452 

6.3 

-.2177 

-.2395 

13.865 

13.5 

5000 

-.2596 

- . 2785 

6.481 

6.3 

-.2104 

- . 2380 

13.878 

13.5 

8000 

- . 2666 

- . 2869 

6.520 

6.3 

- . 2244 

-.2452 

13.895 

13.5 

10000 

- . 2943 

-.3107 

6.648 

6.3 

- . 2565 

- . 2656 

13.940 

13.5 
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Table 4.9: Natural Frequencies of Shafts with Internal Hysteresis 
Damping and Undamped Istropic Bearings 


K 

XX 

= K = 33 

yy 

.14 C = C 

xx yy 

= 0.0 „ H 

= 0.0002 

Spin Speed 
rpm 

Growth 

Factor 

Natural 

Frequency 

Present 

Nelson C123 

(Non -dimensional: sedj 

Present 

Nelson C121 

(From graph and 

Non- d imer,& one ? s ed\ 

Ca) Forward 

Precession 




2000 

.0064 

1 . 0724 

.0018 

1.0621 

6.234 

13.367 

6.434 

13.6 

4000 

.0065 

- 

6.234 

6.434 

1.0738 


13.352 

13.6 

8000 

.0068 

.0018 

6.234 

6.434 

1 . 0754 

1 .0643 

13.351 

13.6 

10000 

.0071 

— 

6.235 

6.4 

1 .0789 

— 

13.324 

13.6 

12000 

.0074 

.0018 

6.232 

6.35 

1 .0813 

1.0665 

13.323 

13.5 

Cb> Backward 

Precession 



2000 

-.2353 
-1 .0526 

-.2546 

-1.0621 

6.465 

13.621 

6.434 

13.6 

4000 

-.2350 
-1 .0564 

- 

6.435 

13.622 

6.434 

13.6 

8000 

-.2324 

-.2563 

6.436 

6.434 

-1 .0602 

-1 .0752 

13.636 

13.6 

10000 

-.2322 
-1 .0832 

- 

6.442 

13.646 

6.45 

13.65 

12000 

-.2426 
-1 .0946 

-.2643 

-1.0852 

6.446 

13.656 

6.50 

13.7 
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Case S: Hysteresis damping, isotropic damped bearings 

The system studied is Case 4 with isotropic bearing dampings 

incorporated as C = C = 0.269. The results are shown in Table 

a a yy 

10. Both the modes are now stable for the entire speed range 
studied. The observation match with those of 112,143. 

In the study of above shaft with internal damping, common 
observations made are listed below: 

(a) The natural frequencies in forward mode of precesion were 
found to be lower and in back ward precesion higher than those 
observed by C 1 2 , 143. This is due to the effect of shear 
deformation which is incorporated in the present work. 

(b ) Forward and backward precesion frequencies were found to be 
very close to each other except when bearing properties were 
anisotropic, in which case they are wide apart. 

(c ) Internal viscous damping is seen to have destabilising effect 
beyond the first critical speed, where as internal hysteresis 
damping destabilises all the forward modes for the entire speed 
range studied. 

(d) Stability of the system is improved by providing damped 
bearing supports. Anisotropic bearing supports also was found to 
increase stability. 

(e) Bearing damping was found to lower the natural frequencies. 


as expected. 
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Table 4.10: 

Natural Frequencies of Shafts 

with Internal 


Hysteresis 

Dampingand Damped 

Istropic 

Bearings 

K 

XX 

= K =33 

yy 

.14 C = C 

xx yy 

0.269 r, H 

= 0.0002 

Spin Speed 
rpm 

Growth 

Factor 

Natural 

Frequency 

Present 

Nelson C123 

(Non-dimensionalised) 

Present 

Nelson 1121 

(Frorr. graph and 
Mon-dimensi'onalised) 

Ca3 Forward 

Precession 




2000 

-.1004 

-.1006 

6.236 

6.35 

-.6821 

-.6838 

13.331 

13.6 

4000 

-.1003 

— 

6.235 

6.35 

- . 6820 

— 

13.330 

13.6 

5000 

-.1002 

-.6817 

-.1003 
- . 6836 

6.234 

13.288 

6.35 

13.5 

8000 

-.1006 

-.6813 

- 

6.223 

13.247 

6.30 

13.4 

9000 

-.1012 

-.1015 

6.212 

6.25 

-.6803 

-.6831 

13.214 

13.35 

Cbi Backward Precession 



2000 

-.1005 
- . 6835 

-.1006 

-.6836 

6.351 

13.624 

6.35 

13.6 

4000 

-.1006 

-.6837 

— 

6.352 

13.625 

6.35 

13.6 

5000 

-.1007 
- . 6838 

-.1013 

-.6843 

6.352 

13.626 

6.35 

13.6 

8000 

-.1007 
- . 6839 

- 

6.352 

13.664 

6.35 

13.65 

9000 

-. 1007 
-.6839 

-.1008 

-.6845 

6.356 

13.697 

6.35 

13.65 
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4.5 SHAFT SUPPORTED IN ELASTIC HYDRODYNAMIC BEARINGS WITH CROSS 
COUPLING 


A steel shaft with following properties supported at both 
ends in similar hydrodynamic bearings is studied for natural 
frequencies over a range of shaft spin speeds. This shaft was 
also studied by C15,173. 


Shaft density 
Young's modulus 
Shaft diameter 
Shaft length 
Bearing L/D ratio 


7833 kg/m' J 
206.8 GPa 
101.6 mm 
1270.0 mm 
0.25, 0.5 


For bearing with L/D ratio as 0.5, the bearing properties 
were obtained from the polynomial functions given by C53. When 
the L/D ratio was 0.25, direct values were not available and were 
extrapolated from graphs for L/D ratios of 0.5 and 1.0 given by 
C53 . 

Results for bearing with L/D ratio of 0.5 are given in Table 
11 and compared with those given by C 15 3 , obtained by FEM. The 
present frequencies are slightly lower than those of C153 as shear 
deformations have been accounted for in the present work. 

Results for bearings with L/D ratio of 0.25 are given m 
Table 12 and compared with those given by C173, using transfer 
matrix method with lumped mass approach. Present frequencies are 
higher than those of 11173. This is expected to a large extent as 
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Table 4.11: Natural Frequencies of Shaft Supported in Elastic 
Hydrodynamic Bearing (L/D = 0.5) 


Bearing Properties: Length = 50.8 mm 

Diameter = 101.6 mm 

Radial Clearance = .0508 mm 

Viscosity of Oil = 6.9 Centipoise 


Spin Speed 
rpm 

Growth 

Factor 

Natural 

Frequency 

Present 

Glassgow C151 

(Mon-dimensi'onali'sed ) 

Present 

Glassgow C 153 

C^on- dimensional! $&d) 

•mm jmL jm. 

-1.4343 

—1 -6809 

1.3208 

1.5508 

2.000 

-1 .6384 

- 

1 . 1684 

- 


-2.2432 

- 

9.6856 

— 


- .8643 

— 

9.4324 

— 


-1.9325 

-2.1215 

3.2623 

3.5642 

4000 

-2.3264 

- 

3.0646 

- 


-2.6864 

- 

9.7356 

— 


- .9852 

— 

9.5358 



- .9367 

-1.0590 

5.0175 

5.2809 

6000 

-2.4643 

- 

4.8964 

— 


-2.9384 

- 

9.8863 

*— 


-1.6987 

- 

9.6634 

— 


- .1678 

- .2398 

5.5639 

5.8406 

8000 

-2.8324 

- 

5.0634 

— 


-3.6047 

- 

9.9264 



-1.9964 

— 

9.8643 



.2736 

.2632 

6.0043 

6.1243 

10000 

-2.6364 

- 

5.6746 

— 


-3.8243 

- 

9.9986 



-2.3246 

- 

9.8885 



.6829 

.6553 

6.2059 

6.4235 

12000 

-2.6046 

— 

5 . 7623 

“* 


-4.0246 

— 

10.6432 



-2.5648 

- 

9.9893 




Table 4.12: Natural Frequencies of Shaft Supported in 
Elastic Hydrodynamic Bearing (L/D = 0.25) 

Bearing Properties: Length = 25.4 mm 

Diameter = 101.6 mm 

Radial Clearance = .0508 mm 

Viscosity of Oil = 6.9 Centipoise 


Spin Speed 

Growth 

Factor 

Natural 

Frequency 

rpm 

Present 

Bansal L 171 

(No; 'dimensionaiised j 

Present 

Bansal C173 

(Xbti-dimensionsiisfid.) 

2000 

-2.2855 
-2.0156 
-3.0343 
- .6326 

-1.5205 
-1.5105 
-2. 1321 
- .4561 

2.0129 

2.2270 

11.5220 

10.6343 

1.6190 

1.7333 

9.6309 

9.9175 

3000 

-2.3434 
-2. 1032 
-3.3667 
- .9643 

-1.8102 
-1.7308 
-2.4656 
- .6317 

2.3449 

3.0726 

10.7360 

9.9725 

2.3155 

2.5493 

9.4927 

9.8732 

4000 

-2.7436 
-1 .6343 
-3.8563 
-1.0432 

-1.9507 
-1 .7272 
-2.7302 
- .7595 

3.0343 

3.9641 

10.3241 

9.4487 

2.9996 

3.3568 

9.3996 

9.8398 

6000 

-2.8667 

-1.4634 

-4.1123 

-1.4836 

-2.0568 
-1.2242 
-3.5790 
- .9243 

5.4321 

6.0343 

10.4321 

9.9642 

4.3601 

4.8314 

9.4011 

9.8398 

7000 

-2.8734 
-1.4343 
-4.6321 
— 1 ■ 6800 

-2.0484 
-1.2241 
-4.0720 
- .9594 

6.2147 

6.9923 

10.8843 

9.9144 

5.0395 

5.3711 

9.6076 

9.7855 

8000 

-2.6324 
- .9343 
-4.7343 
-1.7243 

-1.9806 

- .4164 
-4.5850 

- .9754 

6.6762 

7.7324 

11.4323 

9.8836 

5.7057 

5.7500 

9.9718 

9.7435 

9000 

-2.2343 
- .2367 
-4.8434 
-1.9643 

-1.8992 
- .0781 
-4.8410 

7.0173 

7.4361 

11.8843 

9.7323 

6.3712 

6 . 0605 
10.3920 





FIG 4-4 Whirl Orbit At 6000 rpm (Case i) 



Amplitude of Response Uxio 



Spin Speed A 


Fig 4-5 Amplitude Vs Spin Speed (Case 2) 
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Fig- 4 *7 Amplitude Vfe Spin Speed (Case 3) 
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•1 *£ 3 -4 *5 -6 *7 *8 *9 *0 *1 *2 1*3 

Spin Speed fl 


Fig 4*9 Amplitude Vs Spin Speed (Case 4) 




1.4-10 Whirl Orbit At 6000 r pm [Case 4) 




Fig. 4*11 Amplitude Spin Speed [Case 5) 



CHAPTER V 


CONCLUSIONS 

The fallowing conclusions emerge from the present work: 

1) Transfer finite element method is easily applicable for study 
of flexural vibrations of rotating shops. Results attained 
with fewer elements are quite accurate. 

2) By using four second order differential equations, all types 
of geometrical and natural boundary conditions, including 
bearing cross coupled properties can easily be incorporated 
exactly. 

3) When shear deformation is considered, results obtained using 
cubic polynomials approximation are very accurate with no ill 
condi tioning . 

4) Neglecting shaft mass can lead large errors. The mass of the 
shaft is taken care of in a routine way in the present work. 

5) Effect of shear deformation was found to lower natural 
frequencies in forward whirl and increase in backward whirl. 
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6) Internal viscous damping begins destabilising effect after a 
certain spin speed is reached. Hysteretic internal damping 
has destabilising effect at all speed. 

7) Damped bearing supports increases system stability. 
Anisotropic bearing are also found to improve stability. 

8) For unbalanced shaft system with damping, there is a phase 
shift of response along axial length. 

9) The whirl orbit for isotropic bearings are circular, while 
for anisotropic bearings they are elliptical, as expected. 



APPENDIX 


DERIVATION OF KINETIC ENERGY 

The shaft cross section has a rotation 8 about x— axis, 8 

x y 

about y-axis, 8 about s-axis, translation U along x-axis and V 
along y-axis as mentioned in section (2.1> - The kinetic energy 
due to rotation is obtained as fallows: 

Consider a frame of reference abc attached to shaft cross 
section such that c-axis is normal to the cross section, a,b axes 
are in the plane of cross section and mutually perpendicular to 
each other and c-axis« Then the rotations of cross sections, 
relative to xys frame are also the rotations of frame abc relative 
to xys- Initially the axes of abc frame are considered to be 
coincident with xys frame. The final position of frame, abc 


relative 

to 

xys is 

determined 

by the following 

sequence of 

rotations 

: 






1 ) 

A 

rotation 

by 8 

J X 

about 

x-axis gives a 0 b 7 c ? 

(Fig. A- 1 a > 

2) 

A 

rotation 

by 8 

y 

about 

b 7 ~axis gives a^b^c^ 

(Fig. A-ib) 

3) 

A 

rotation 

by 6 ^ 

about 

c^-axis gives abc 

(Fig. A— 1c ) 


Referring to Fig. (A-ia), the position of a^b^c^ frame after 
the rotation of 6 relative to xys frame as 
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a 2 



' 1 

0 

0 


f \ 

X 

b 2 

Y 

= 

0 

cos 6 

z 

sin 6 
- z 

1 

y > 

C 2 J 



0 

-sin 6 

z 

cos 6 

z J 


s 


(A— 1 . 1 ) 


after 


Similarly, the position of frame a 1 b 1 c 1 
the rotation 6 is, Fig.(A-lb). 

y 


relative to a^^c^ 


f *v 

a. 



COS 

6 

0 

-sin 

e 1 


f \ 

a_ 

1 




y 



y 


2 

b l 

Y 

5= 

0 


l 

0 



b 2 ' 

c.. 



sin 

e 

0 

cos 

6 


c_ 

1 - J 



L 

y 



y J 


K 2 


(A-1.2) 


Finally, the position of frame abc relative to 
after the rotation & : is. Fig. (A-lc). 



frame a^b^c^ 


( A— 1.3) 


The final position of abc after the above rotations is shown 
in Fig. (A-ld) . The frame abc, has the angular velocities 6^ 
along a^, 6^ along b^ and 6 __ along c. The component of 6^ along 
the directions a , b and c can be obtained by the following 
transformations using eqns. (A-1.2) and (A-1.3) as 
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Similarly the components of 6 ' along the directions a,b and c 
are obtained using eqn. (A-1.3) as 
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(A-1.5) 


The angular velocity vector 8^ is along c direction, thus no 
transformation is required. Hence the angular velocities along a, 
b and c directions, due to the rotations & 8 and 6^ are 


obtained 

using eqns (A-l 

. 4 ) and 

<A— 

1.5) 

as 
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As 

the directions 

a,b 

and 

c 

are 

also 

the 

principle 

directions, the kinetic 

energy 

due 

to 

translation 

and rotation 


becomes , 


Te = J |pA(U 2 +V 2 > + + I p ^ j ds (A-1.7) 


Substituting eqn. (A-1.6) the expression for kinetic energy 


becomes 
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Assuming & . 6 to be small and noting that 6 - O <the 

x y z o 

speed of the shaft), the above expression simplifies to 
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(c) For isotropic bearing properties, the whirl orbit is circular 
and when bearings are anisotropic the whirl orbit is elliptical as 
expected . 

(d) With damping present, there is a phase shift of response 
along axial length. 

(e) The response is symmetric about mid plane for all cases 
except Case 5, when the response at a point to the right of mid 
plane is lower than the corresponding point at left side as 
expected, as bearing at right hand end is stiffer than the left 


end one. 



